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ABSTRACT 

We have monitored the position of the cool Y dwarf WISEPA J182831. 08+265037.8 
using a combination of ground- and space-based telescopes and have determined its 
distance to be 11.2^^ pc. Its absolute H magnitude, M# = 22.21+Q22 mag, suggests a 
mass in the range 0.5-20 Mj up for ages of 0.1-10 Gyr with an effective temperature in the 
range 250-400 K. The broad range in mass is due primarily to the unknown age of the 
object. Since the high tangential velocity of the object, 51±5 km s _1 , is characteristic 
of an old disk population, a plausible age range of 2-4 Gyr leads to a mass range of 3-6 
Mj up based on fits to the (highly uncertain) COND evolutionary models. The range 
in temperature is due to the fact that no single model adequately represents the 1-5 
/im spectral energy distribution (SED) of the source, failing by factors of up to 5 at 
either the short or long wavelength portions of the spectral energy distribution. The 
appearance of this very cold object may be affected by non-equilibrium chemistry or low 
temperature condensates forming clouds, two atmospheric processes that are known to 
be important in brown dwarf atmospheres but have proven difficult to model. Finally, 
we argue that there would have to be a very steep upturn in the number density of 
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late type Y-dwarfs to account for the putative population of objects suggested by by 
recent microlensing observations. Whether WISE 1828+2650 sits at the low mass end 
of the brown dwarf population or is the first example of a large number of "free-floating" 
planets is not yet known. 

Subject headings: brown dwarfs infrared:stars - individual (WISEP J182831. 08+265037.8) 
- Astrometry - Parallaxes - proper motions - solar neighborhood 



1. Introduction 



A long term goal of infrared sky surveys has been the d etection and characterization of sub 



stella r objects of progressively lower temperature and mass (jKirkpatrick, Beichman fc Skrutskie 



19971 ). The Two Micron All sky surv ey (2MASS), DEep Near Infrared Southern Sky Survey (DE - 
NIS) and Sloan Digital Sky Survey dskrutskie et alJbood : Epchtein et aDll997l : Kork et alJ I2OOOI ) 
identified the first bona fide field brown d warfs while coronagraphic studies found the fir st bona 
fide brown dwarf companion, Gliese 229B (jNakajima et al.lll995l : lOppenheimer et al.lll998l ). These 
studies led to the acceptance of the first new spectral classes to be named in over 50 years, the 
L and T dwarfs (IKirkpatrick et al.lll999l ) characterized by effective temperatures, T e ff, as low as 
600 K. Studies of these objects suggested the potential for yet a third class of even colder objects, 
the Y dwarfs, with temperatures between 300-600 K. Many exa mples of L and T d warfs have now 
been discovered and their spectral classification is now a fine art (IKirkpatrick 1120051 ) . But members 
of the Y spectral class were elusive until the launch of the Wide-field Infrared Survey Explorer 
(WISE; Wright et al 2010) w hich identified over 100 brown dwarfs with lo w effective temperature, 



including the first Y dwarfs. (IKirkpatrick et al.ll201ll ; ICushing et al.l 12011 



In parallel with the WISE survey, targeted imaging of stars has led to the discoveries of 
planetary mass obj ects orbitin g nearby stars, yieldi ng; both young planet s, e.g. HR8799 bcd e 



dMarois et al.l boosf l. /? Pic b Lagrange et~al1 \20od ) and 2M1207-3 932b dchauvin et al 



as well as a companion to a nearby white dwarf, WD 0806-661 b (ILuhman et al 



2011 




Thermal emission has also been detected from Jovian mass planets transiting their parent stars 
(IDeming et al.l 12006k ICharbonneau et al.l 120051 ) . Finally, microlensing studies may have identified 
a population of Jovian-mass objects with no obvious association to any star and which th us must 



either be located tens of AU from their hosts or be free floating planets (jSumi et al.l 12011 



Essential to understanding the nature, physical properties and evolutionary status of these 
objects is the determination of such properties such as mass, age, composition, and bolometric 
luminosity. This article describes the first results of a program to use direct imaging to determine 
distances to these ultra-cool brown dwarfs using parallax measurements. With secure distance 
measurements it will be possible to estimate the local space density of Y dwarfs and to understand 
better their physical properties. 
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We have taken a sample of twenty objects from the entire all-sky WISE survey that are 
confirmed to have spectral types (or types estimated from WISE photometry) of T8 and later and 
have initiated a program of follow-up photometry from ground- and space-based teles copes. From a 



compa r ison of apparen t and inferred absolute magnitudes for objects of low T e ff , e.g. iBaraffe et al 



(|2003l ); iBurrows et al.l (j2003l ), we estimate that these objects are within 20 pc of the Sun and so 
shoud have large astrometric signatures. We are using the Hubble Space Telescope (HST), the 
Spitzer Space Telescope, and the Keck-II telescope to obtain images with the goal of determining 
parallactic distances. This paper focuses on the first result of this program, the determination of 
the distance to WISEPA J182831. 08+265037.8, hereafter, WISE 1828+2650 which was identified 
as a poten tial late-type brown dw arf on the basis of its extreme color, W1-W2 ([3.5/im]-[4.6/im]) 
> 4 mag (IKirkpatrick et al.l l201ll ). S ubsequent near-infra red observations from HST led to its 
classification as the prototype Y dwarf (ICushing et al.ll201ll ). and fin ally, on the basis of a br oader 
sample of Y dwarf discoveries have revised the classification to >Y2 (IKirkpatrick et al.ll2012l ). 



2. Observations 

The determination of high accuracy positions for cool, low luminosity brown dwarfs is extremely 
challenging since they are faint for 8-10 m ground-based telescopes and HST in the near-infrared 
and bright in the mid-IR where they can only be observed with space telescopes of small aperture 
(0.4-0.85 m). As a result we have had to piece together positional information with a variety of 
telescopes over a range of wavelengths using individual measurements (Table [1]) with uncertainties 
in the range of 10-50 mas due to either low signal to noise ratio or large beam sizes. A parallax 
determination will be possible only if these objects prove to quite close-by. Fortunately, this appears 
to be the case for at least WISE 1828+2650 where we are able to determine a parallax of ^100 
mas. With such a large astrometric signature, relatively coarse positions suffice and many of the 
instrumental effects that must be taken into account for more distant objects can be ignored. 
Table [T] lists the observations from WISE and other observatories. Photometric data are presented 
in Table [3l For WISE we report an average of the 4.6 /im magnitudes from the two epochs and 3a 
upper limits in the other bands. For comp arison with models we have converted the magnitudes 



to flux densities using the zero points from lWright et al.l (120101 ). but because of the unknown and 



extremely non-blackbody-like nature of brown dwarf spectral energy distributions (SED), we have 
not color-corrected these flux densities. 



2.1. HST Observations 

WISE 1828+2650 was observed with HST in the WFC3-IR F140W filter in coni unction with a 
grism spectrum which was the basis of this object's identification as a >Y1 dwarf by lCushing et al 



(1201 ll ) (Figure [T]). The data consist of two pairs of individual 78 sec exposures taken at positions 



offset by (Aa, AS) =0.82" and 5. 67" to reduce the effects of cosmic rays and the undersampling of 
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the individual frames. The pipelin e corrects for the geometric d istortion of the WFC3-IR camera 



to a level estimated to be ~ 5 mas (Kozhurina-Platais et al 



20091 ) which is of the same order or less 



than the extraction uncertainties of the faint target. Over 200 sources were extracted from the Space 
Telescope Science Institute's (STScI) "MultiDrizzle" mosaic using the Gaussian-fitting IDL FIND 
routine to determine centroid positions and the APER routing with a 3 pixel radius for photometric 
measurements. The Full Width at Half Maximum (FWHM) fo r these undersampled data is 0.26" 
consistent with STScI analyses fjKozhurina-Platais et al.ll2009l ). The data were calibrate d using 



the z ero-point for Vega magnitudes, ZP(Vega)=25.1845 magj, from the WFC3 Handbook (IRajan 



2010). A strong source denoted Source A appears to be slightly extended and is located roughly 



5" W and 3" S of WISE 1828+2650. This object is seen in the Keck, HST, and Spitzer images 
and will be used as a local astrometric reference in subsequent analysis. As discussed below, the 
extracted position of Source A is consistent between the Keck and HST measurements to within 
10 mas which is small compared to the other uncertainties. 



2.2. Keck NIRC2 Observations 



WISE 1828+2650 was observed in the H-band using the Keck-2 telescope (July 2010) using 
NIRC 2 with the laser guide star Adaptive Optics (AO) system flWizinowich et al ESI Ivan Dam et al 
20061 ) and a tip-tilt star located 50" away. The wide- field camera (40 mas/pixel scale; 40" field 
of view) was used to ensure an adequate number of reference stars for astrometry. A sequence of 
images with a total integration time of 10 80 sec was obtained at airmass of 1.0-1.1 with 2 // -10 // 
offsets in right ascension and declination (IKirkpatrick et al.l 1201 ll ) . The individual images were 
sky-subtracted with a sky frame created by the median of the science frames and flat-fielded with a 
dome flat using standard and custom IDL routines. Individual images were "de-warped" to account 
for optical distortion in the NIRC2 camera as described below. The reduced images were shifted to 
align stars on a common, larger grid and the median average of overlapping pixels was computed 
to make the final mosai c. The target was detected at a relatively low SNR of ^5. As reported in 
Kirkpatrick et al.l (120121 ). the source was also weakly detected in the J band with NIRC2 at this 
epoch (Table E]). 

The source was re-observed in 2011 (Figure [2]) and twice in 2012 following the same procedure 
but for a longer total integration times (3960, 2880 and 2880 sec, respectively), yielding SNR^ 10 
and with smaller offsets (2-4") to reduce the effects of optical distortion on the final image. The 
airmass for the Oct 2011 observations ranged from 1.1 to 1.4 with a mean value of 1.25. In 2012 
July the airmass values ranged between 1.04 to 1.15 with a mean value of 1.09. In 2012 October, 
the seeing and photometric conditions were relatively poor with an average airmass of 1.2. The 



1 A\\ the photometric measurements reported herein were made using this routine from the IDL ASTRO library, 
http:/ /idlastro. gsfc.nasa .gov/contents. html] A number of other IDL routines are taken from this library as well. 

2 http: / / www.stsci.edu/hst / wfc3 / phot_zp_lbn 
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typical Point Spread Function (PSF) in the final mosaicked images from the first three epochs varies 
slightly across the field, but at the center of the field close to WISE 1828+2650, the Full Width 
Half Maximum (FWHM) of the AO-corrected core is ^3-4 pixels, or ^120-160 mas. In October 
2012, the image quality was about a factor of two worse as is reflected in the poorer astrometric 
solution described below. Sources were extracted using the IDL APER routine with an aperture 
radius of 5 pixels. To calibrate these data photometrically we used obs ervations of the W ISE 
1828+2650 field obtained with the Wide Field Infrar ed Camera (WIRC) (Iwilson etalJ booah at 
the Palomar 5-m telescope with a MKO-NIR H filter f|Tokunaga et al.ll2002l ). The 8.7' field of this 
camera was sufficient to detect ^350 2MASS stars surrounding WISE 1828+2650 as well as 6 stars 
seen in the Keck field. While the Palomar image was not deep enough to detect WISE1828+2650 
itself, we were able to calibrate the Keck field using the calibrated W IRC magnitudes of the 6 stars 
in common. Combined with the 2MASS calibration uncertainties (jCarpenter I l200ll ). the overall 
calibration uncertainty is ^5% which is small compared to the uncertainty in the average magnitude 
for the WISE target itself (~ 9%). There is no evidence for variability between the different epochs. 
No photometry was attempted for 2012 October. 

The source positions obtained from the Keck images were corrected for the effects of differ- 
ential refraction relative to the center of the field using meteorological conditions available at the 
CFHT telescope weather archival to determine the index of refraction co rrected for wavelength, 
l ocal temper ature, atmospheric pressure and relative humidity (jLang Ill983l ) and standard formulae 
(jstonel[r996l ). As discussed in Appendix I, for the small field of view of the NIRC2 images and the 
relatively low airmasses under consideration here (<1.5), the first order differential corrections we 
do apply to the source positions are small, <10 mas across the ±20" field and proportionately less 
at smaller separations. 

The effects of optical distortion in the wide- field NIRC2 camera were corrected using a dis- 
tortion map derived by comparing Keck data of the globular cluster M15 taken immediately after 
the WISE 1828+2650 observations in 2 011 and July 2012 to HST/WFC3 observations of M15 
fjAnderson et al.ll2008l : lYelda et al.ll2010l ). Details of this distortion mapping are described in Ap- 
pendix I but the correction amounts to < 1 pixel (40 mas) across most of the array and up to 
2 pixels at the edges of the array. At the array center where the WISE 1828+2650 observations 
were made, the maximum uncorrected distortion for separations up to 5" is < 20 mas. After our 
correction procedure the residual distortion errors are less than 10 mas over the entire field. As 
a test of our ability to co-register data from individual Keck images, we examined the positional 
dispersions of 15 sources in each apparition with respect to their average positions. For the first 
three Keck epochs the average dispersion in right ascension and declination is ±6 mas. The last 
epoch was not considered in this analysis due to poor observing conditions. The ability to register 
these images to an external reference, i.e. the sources in the HST image, is discussed below (§3]). 



3 htt p : / / mk wc . ifa. hawaii . edu / archive / wx / cfht / 
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2.3. Spitzer Observations 

The Spitzer Space Telescope observed WISE 1828+2650 using Director's Discret ionary Time 



and at three subsequent General Observer epochs. The IRAC camera fjFazio et al.ll2004l ) was used in 
its full array mode to make observations at 3.5 (Channel-1) and 4.5 /im (Channel-2) in the first epoch 
and only at 4.5 /im in the subsequent epochs. We analyzed post-BCD mosaics from the Spitzer 
Science Center (SSC) to make photometric and astrometric measurements, extracting sources using 
a 3 pixel radius aperture and normalized using SSC-recommended aperture correction^. WISE 
1828+2650 was detected at both 3.5 /im (SNR~ 10) and at 4.5 /im (SNR>40; Table EJ). Source A 
is seen as a weak but definite detection with about one-sixth the brightness of WISE 1828+2650 
(SNR=9). We fitted the pixel data for these two objects simultaneously using Gaussian profiles to 
determine accurately the source centroids. The fitting uncertainties for the two objects in the 4 
Spitzer observations are 10-20 mas for WISE 1828+2650 and 120-130 mas for Source A (Table ED- 
There was no evidence for variability in the brightness of WISE 1820+2650 between the four epochs 
at the <6% level. Since Source A appears extended at shorter wavelengths, there is a concern that 
position of Source A might shift with wavelength. There is no evidence for this effect in our data. 
The Spitzer position for Source A is self-consistent within 15 mas over 4 epochs and consistent with 
the Keck and HST positions to within (Aa, A5)=(60,10) mas which is small relative to the fitting 
uncertainties for Source A of 120 mas. 

The IRAC camera has significant optical distortions over its entire field that are corrected in 
the SSC pipeline to an average residual value of 100 mas over the whole flelc0. We have used 
-20 HST stars observed within ±90" of WISE 1828+2650 as probes of the residual distortion 
as described above for the Keck observations. As a test of our ability to co-register data from 
individual Spitzer images, we examined the positional dispersions of 22 sources in each apparition 
with respect to their average positions. The average dispersion in the four Spitzer epochs is ±60 
mas in right ascension and declination. The ability to register these images to an external reference, 
i.e. the sources in the HST image, is discussed below (^3j). These values are less than the quoted 
100 mas uncertainty from the SSC, in part because we are confining our observations to the central 
part of the IRAC array. The total uncertainties in Table [5] reflect the systematic uncertainties 
combined with a brightness-dependent term for WISE1828. In the case where we consider only 
offsets from WISE 1828+2650 to Source A, the distortion effects will be considerably smaller given 
the ^10 pixel distance between the two objects. From examination of the polynomial expression to 
the distortion correction given by the SSC, we estimate the distortion component of the uncertainty 
in the separation between the two sources to be <10 mas, negligible compared to the uncertainties 
due to the source fitting (Table [6]) . 



4 http:/ /irsa.ipac. caltech.edu/data/SPITZER/docs/irac 

5 http: / / irsa.ipac.caltech.edu / data/SPITZER / docs / irac / iracinstrumenthandbook /26/ 
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3. Astrometric Data Reduction 



Following source extraction, we put observations from all telescopes onto a common astro- 
metric frame. This was accomplished in a two-step process. First, to ensure consistency with 
the WISE positions which are based on the 2MASS astrometric frame (j Wright et al.l bold ), the 
HST positions were refined using 11 2MASS sources to solve for the RA and DEC zero points 
and image rotation. We adopted the no minal STScI pipeline plate sca le of 0.12826 arcsec/pixel 
which is known to an accuracy ~ 0.01% (jKozhurina-Platais et al.ll2009l ). This value agrees with 
the plate scale derived from the 2MASS sources to within <0.02%, but we adopted the better 
established pipeline value. Saturation effects on the astrometric centroid for the faint 2MASS 
sources (H>15 mag) are small (<< the 0.128" HST pixel) compared to the overall 80 mas uncer- 
tainty in the absolute 2MASS reference frame (ISkrutskie et al.ll2006l ). Additional evidence that 
the saturation effects are relatively small comes from the fact that our HST F140W photometry 
recovers the 2MASS brightness of these 11 objects to within an average ±0.2 mag where we com- 
pared the HST data with the simple average of the 2MASS J and H magnitudes. The fit to the 
HST/2MASS sources yielded a (Aa , AS) (initial HST-2MASS)=(-0.50±0.02,0.15±0.05)" shift rela- 
tive to the initial pipeline solution where we adopted the positional uncertainties from the 2MASS 
catalog to weight the individual sources (70-350 mas with an average value of 180 mas) in the 
fitting procedure. We believe the 2MASS positional zero points to be correct given the agreement 
between 9 sources within 2! of WISE1828+2650 appearing in both the 2MASS and the UCAC3 
catalogs (jZacharias et al.ll2010l ). The average positional differences between these two catalogs is 
(Aa, A5)(UCAC3-2MASS)=(0.12±0.01,0.03±0.02) // . 

Sources seen in the HST image and seen in either the Keck, Spitzer, or 2MASS data are 
presented in Table [2 The positional errors associated with th e HST data come f rom a combination 
of a brightness-dependent term, a S NR = FWHM/(2 x SNR) (jMonet et alJboid ) where SNR is the 
signal-to-noise ratio of a particular measurement, and two terms term depending on distance of the 
source, R, from the center of the HST field, a rot = Rd6 and a sca i e = i?d(scale) /scale. The latter 
two terms account for uncertainty in the rotation angle and plate scale. The rotational uncertainty 
is d6 = 0.028° = 4.8 x 10 -4 rad (Tabled]). For sources at the edge of the field the positional 
uncertainty can be as large as 90 mas, but within 20" of the center of the Keck field, for example, 
the typical uncertainty of an HST source is ^5 mas. 

Parameters of the 2MASS/HST fit are given in Tabled! As shown in the first entry in Table [11 
the standard error (cr/y/~N) in this local frame is approximately 4 mas in each axis. To this local 
registratio n uncertainty must be added the overall astrometric accuracy of the 2MASS survey of 
^80 mas (jSkrutskie et al.ll2006l ). However for the differential measurements of interest for the 
parallax determination, this uncertainty is not important. The new solution was applied to all the 
HST sources, of which 13-22 appear in various of the Keck and Spitzer images (Table ED- These 
refined HST positions were subsequently used to set the astrometric reference frame for all other 
images. 
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We registered the four epochs of Keck imaging to the 2MASS / HST astrometric grid using 
between 13-21 stars, solving simultaneously for field center coordinates, rotation, and pixel scale. 
The solutions yielded offsets relative to the HST frame of <2.5 mas, with a single-axis standard 
error ^3 mas, and fitting errors of < 0.003° and 0.02% in field rotation and plate scale, respectively. 
These results are summarized in Table HI We averaged the Keck positions of the HST sources for the 
first three epochs and compared these to the input HST positions, giving (Aa , A5)(Keck-HST) = (- 
9+5, 5+8 mas) with a dispersion of 24 mas predominantly independent of source brightness. We 
take this dispersion to be representative of the residual uncertainties in the Adaptive Optics PSF 
across the field, uncertainties in the rotation and plate csale of the HST data (Table EJ), residual 
optical distortion in both the Keck and HST data, and residual refraction in the Keck data. To 
establish the uncertainty for a particular object we combine this observed dispersion in a root-mean- 
square sense with the brightness-dependent term, cjsnr- This latter term is typically small, <10 
mas. As noted above, the 2012 October Keck data were obtained under poor conditions resulting 
in a dispersion of 50 mas with respect to the HST sources. 

For Spitzer astrometry we used the 4.6 fim (Band 2) Post-BCD mosaics because the source is 
brighter at this wavelength compared to 3.5 fim and also because of the more symmetric and better 
sampled PSF at the longer wavelength. We matched HST positions of ^20 stars (depending on the 
intersection of the two fields of view), solving for RA and DEC zero points, rotation angle, and pixel 
scale. As expected, the Spitzer scale proved to be highly stable over the four epochs of observation 
and to differ from the SSC value by <0.05%. We adopted the nominal plate scale for our analysis. 
We determined the rotation angle to < 0.02° and the positional zero point to 10 ^20 mas relative 
to the HST frame. We averaged the Spitzer positions of the HST sources for the 4 epochs and 
compared these to the input HST positions, giving (Aa , A5)(Spitzer-HST) = (7±9,-5±ll) mas with 
a dispersion of 50 mas predominantly independent of source brightness. We take this value to be 
representative of the residual uncertainties in PSF and distortion in both the Spitzer and HST data 
that must be applied to all sources. We adopt the slightly larger repeatability measure of ±60 mas 
obtained in the previous section as the astrometric uncertainty in the Spitzer data; the SNR-related 
term is negligible by comparison. 



Results 



4.1. Determination of Parallax and Proper Motion 



Table E lists the positions of WISE 1828+2650 on the J3ST/2MASS astrometric grid while 
Table [6] lists the offsets of WISE 1828+2650 relative to source A. Both sets of right ascension and 
declination d ata were fitted to a model incorporating proper motion and parallax (jSmart 1119771 : 



Green 111983 ): 
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ol =ao + n a (t — Tq)/cos(5 / ) 
S f = 5 + w(* - T ) 
a(t) — a + 7T^X{t)sina — y(£)cosc/^ /cos 5' 

5{t) =5' + tt( X(t)cos a 1 sin 5 r + Y(t)sin a' sin 




(i) 



(2) 



where («o,5o) are the source position for equinox and epoch To =J2000.0, ji a j are proper motion 
in the two coordinates in "/yr, and tt is the annual parallax. The coefficients X(t), Y(t), and 
Z(t) are the rectangular coordinates of a terrestrial or earth orbiting telescope (WISE, Keck, HST) 
or Spitzer in its Earth-trailing orbit as seen from the Sun. Values of X, Y, Z for the terrestrial 
or Earth-orbiting observatories are taken from the IDL ASTRO routine XYZ. Spitzer values of 
X, Y, Z are obtained from the image headers provided by the SSC. Equations (1) and (2) are 
solved simultaneously using the Mathematica routine NonLinearModelFit incorporating appropriate 
uncertainties for each datapoint. Figure [3] shows the fit of the data to this model. Figure [H shows 
the fit to parallactic ellipse with the effects of proper motion removed. 

We developed two solutions, one using best estimates of the time-varying positions of WISE 
1828+2650 in the J2ST/2MASS astrometric frame (Method-1; Table Ej) and the other using offsets 
between the source and the nearby object denoted Source A (Method-2; Table Ej). The Method-2 
position offsets use only the rotation angle and plate scale derived from the full astrometric solution, 
and should thus be relatively immune to residual distortion and refraction effects that apply across 
the larger field. On the other hand, the only Method-2 observable is the offset which is compromised 
by Source-A being both slightly extended in all the images and blended with WISE 1828+2650 in 
the Spitzer images. 

The two solutions give consistent results (Table Ej), yielding parallaxes of 103+16 mas and 
79+12 mas and corresponding to distances of 9.7t.\'i pc (Method-1) and 12. 6^9 pc (Method-2). 
To generate our final kinematic estimates we combined and fitted simultaneously the Method-1 and 
Method-2 datasets, Method-(l+2). While the input imaging data are the same in both cases, the 
sources of uncertainty are different and furthermore, the blending in the Method-2 offsets lowers 
the significance of the Spitzer results and thereby halves the number of epochs with precise mea- 
surements. Figure [5] shows histograms for the parallax and distance solutions from Method-(l+2) 
resulting from a 10,000-iteration Monte Carlo calculation wherein the source positions were varied 
according to their Gaussian uncertainties and solutions recorded at each step. While the parallax 
distribution is consistent with the nominal average value and its associated uncertainties, the dis- 
tance distribution shows a tail to greater distances expected due to its inverse relationship to the 
parallax. We derived the asymmetric la distance uncertainties from the Cumulative Distribution 
Function of the errors (Figure |5b)- Our best estimate of the distance to WISE1828+2650 is 11. 2^^ 
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pc with a x 2 — 36.3 with 31 degrees of freedom. A model of proper motion plus a non-zero parallax 
provides a much better fit to the data than a proper motion-only model which yields x 2 — 135 with 
32 degrees of freedom for the Method-(l+2) dataset. On the other hand, there is no evidence that 
a more complex model, e.g. one including the presence of a perturbing second body in the system, 
is required. 

Finally, we note that the parallax determination here is a relative one with respect to nearby ref- 
erences stars (not galaxies, except for Source A). Thus the motions of the reference fr ame can intro- 



duce a bias into the parallax determination. Since this bias is at the level of 0.5-2 mas (jDupuy fc Liu 



20121 ). we can neglect this correction relative to the 100 mas signal of WISE 1828+2650. 



4.2. Spectral Energy Distribution 



WISE 1828+2650 is one of the coolest (latest) Y dwarfs yet identified. Available photometry 
of WISE 1828+2650 covers the 1-5 fim range and can compared against the predictions of brown 
dwarf atmospheric and evolutionary models (IBaraffe et al.l 120031 ). Figure [6] shows a number of 
models fitted to the observations as described below. Unfortunately, none of these models provide 
a good fit to the full range of the data. The COND models pass through either the two short 
wavelength points or the 3-5 fim points, but not both. The over- or under-estimates are significant, 
up to factors of 5, which suggests that the models have omitted one or more important physical 
processes in th e atmosphere, e.g. n on-equilibrium chemistry or new absorbers, that dramatically 



affect the SED. iMarley et al.l (120101 ) highlight the importance of clouds in the L/T dwarf transition 
which greatly change the near-IR (J, H and K band) appearance of these warmer objects. The 



onset of still lower temperature condensates, e.g. Na2S producing s ulfide clouds (|Morlev et al 



20121 ). may play a role in the appearance of the Y dwarfs (see however iBurrows et al 



(12003^ 



While acknowledging the failings of the models in this uncharted range of effective temperature, 
we will use them to make preliminary estimates of the properties of WISE 1828+2650. We take 
advantage of our knowledge of the distance to the source and use absolute magnitudes rather than 
relative fluxes to constrain further the evolutionary status of the object. 



Magnitudes corresponding to the dust-free COND models (IBaraffe et al.ll2003l ) were calculated 
using bandpasses matched to HST, WISE, Spitzer, and ground-based filters. Models were first 
selected to match either the absolute H magnitude M# = 22.21^Q22 ma § or M[4.5] = M.42^Q2o 
mag where the uncertainty in the absolute magnitude is dominated by the distance uncertainty. 
Selecting models matching either absolute magnitude led to two classes of object. Matching the 
H band absolute magnitude leads to atmospheres with a narrow range of effective temperatures 
(T e ff = 275 ± 40K) but with a broad variation of ages and masses from (0.5 Mj up , 0.05 Gyr) to (8 
Mj wp , 10 Gyr). Matching the [4.5/im] absolute magnitude leads to hotter, more massive objects 
(T e ff = 450 ± 40K) with a range of ages and masses from (2 Mj up , 0.1 Gyr) to (30 Mj up , 10 Gyr). 
Matching the absolute magnitudes in F140W, H, and the two Spitzer bands simultaneously, each 
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with an uncertainty of ^0.25 mag, yields a best fitting model with (0.1 Mj up , 0.1 Gyr), T e ff = 285K, 
and a x 2 — 95 for 2 degrees of freedom (4 data points and two parameters). The variation of the 
quality of fit as a function of mass and age (Figure [Tj) shows a minimum for a young, low mass 
object, but a broad valley tending to older, more massive objects. We emphasize the poor quality 
of the model fits by pointing out that the contour values in Figure [7] show the logarithm of the x 2 
value. 

Finally, we examined the pos sibility that the SEP could be the result of a combination of two 
coeval bodies in a binary system fjLeggett et al.ll2012l ). Fitting the 4 absolute magnitudes yields a 
pair of equal mass objects with (0.5 Mj upi 0.05 Gyr), T e ff = 285K, and x 2 — 55. However, since 
none of these models provides a particularly good fit to the data across the entire 1-5 fim range, the 
estimated physical parameters should be taken with a grain of salt. What all the solutions have in 
common is a low T e ff ~ 250 — 450K and a mass range of a 0.5-20 Mj up depending on the unknown 
age of the system ( §5.2j) . 

Other models lead to a similar range of predicted properties for WISE 1828+26 50 depending 



on the wav elength used for the comparison and the assumed mass and/or age. The Marley et al 



2010) and 



Burrows et al 



models also suggest an effective temperature ~ 300-400 K which 



can be associated with a broad range of masses a nd ages with no sin g le mo del providing a good 



simultaneous fit to all of the 1-5 jim observations. ISpiegel fc Burrows (120121 ) calculate models for 



both "hot start" and "cold start" formation scenarios ( Marley et al 



20071 ) for ages less than 0.1 



Gyr. Matching our data to the predicted absolute H band models for either the "cold start" or 
"hot start" cases yields 1 Mj up objects with ages of 0.02-0.03 Gyr. Matching the predicted M 
band models (only an approximation to our Spitzer filter) yields objects in the range (1-5 Mj up , 
0.007-0.1 Gyr). Differences between the COND models and the Spiegel & Burrows models are 
relatively minor in the area of age overlap. 



5. Discussion 

5.1. Physical Properties of WISE 1828+2650 

The nature of WISE 1828+2650 is ambiguous for two reasons. First, the failure of the models 
to reproduce successfully the SED and absolute magnitudes makes it hard to draw any quantitative 
conclusions. Second, even with improved models, e.g. ones including clouds or new low-temperature 
absorbers, there would remain a family of models with roughly similar T e ff that are degenerate in 
mass and age. 

WISE 1828+2650 is simila r in appearance to t he brown dwarf companion to the ~ 2 Gyr 



old white dwarf WD 0806-661 fjLuhman et al.l 120121 ) located 19.2 pc away. The twofold greater 



distance to this system has so far prevented a detection of WD 0806-661 B in the 1-2 /xm region. 
As a result, we do not yet have a comparable SED, but the near- to mid-IR colors are comparably 
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red with J-[4.5]=9.25±0.35 (IKirkpatrick et all 12012 ) for the WISE object and J-[4.5]>7 for WD 
0806-661 B. The 3-5 fim colors are also similar, [3.6]- [4.5]= 2.77+0.15 mag for WD 0806-661 B and 
[3.6]- [4.5]= 2.81+0.06 mag for WISE 1828+2650, although neither of these objects are as red as 
other Y dwarfs. With an age estimate for the white dwarf i t is possible to narrow the mass range 
for its companion based on its colors. iLuhman et al.1 (|2012l ) suggest a mass of 6-9 Mj up which is 
consistent with our estimates for WISE 1828+2650 in that age range. 

An important peculiarity of WISE 1828+2650 is its bright absolute magnitude, Mw2 = 
14.15^o 21 ma S f° r its presumed Y2 spectral type. Compared with WD 0806-661 B with an abso- 
lute magnitude of M[ 4 5 ]=15.5 mag, WISE 1828+2650 is > 1 mag more luminous. More troubling, 
however, is the fact that compared with other Y dwarfs with reasonably secure parallaxes, WISE 
1828+2650 is > 2 mag brighter than other Y dwarfs in Mw2 and, to a lesser extent, brighter than 
some of the earlier Y dwarfs in Mh • Figure El which illustrates this problem, is an upda ted version 
of the absolute magnitude vs. spectral type figures ap pearing in IKirkpatrick et al.l (120121). This new 



Tinnev et all (12012 ). 



figure uses recently published Y dwarf parallaxes by iMarsh et al.l (120 12l ) and 
The locations of WISE 1828+2650 and WD 0806-661 B are also shown. The spectral type of WD 
0806-661 B has been estimated from its effective temperature since no spectrum is available and 
only an upper limit is available for the H-band magnitude. However, in both cases it is apparent 
that the two objects are similarly deviant from the trend set by earlier Y objects. 

A simple explanation of at least part of the ^1.5 mag difference in absolute M[ 4 5 ] between 
WISE 1828+2650 and WD 0806-661 B is that WISE 1828+2650 is an unresolved, equal mass 
binary which would increase the total system brightness by 0.7 mag, close to the difference with 
WD 0806-661B. Any such companion to WISE 1828+2650 would have to have an orbit < 0.5 AU 
(<0.05 // ) to be consistent with the HST and Keck data and if present would eventually manifest 
itself in either the astrometric data or future spectroscopic data. No such evidence for binarity 
yet exists. In any event, such an explanation does not account for the significant deviation from 
the trends established for earlier spectral types. Some other explanation must be sought for the 
differences in absolute magnitude at 3-4 fim. 



5.2. Age of WISE 1828+2650 

We cannot eliminate solutions corresponding to young, low mass objects for WISE 1828+2650. 



are at le ast three M stars y ounger than 
(6.3 



pc; 



c; Kalas et al. ( 


2004 


)), GJ 896 B 


(jllenry et al. lb 


!006l: Riedel et al. 



Riedel et al.l (120111 )). and most recently AP Col (4.6 pc) (IHenry et al. 
20111 ). A number of brown dwarfs have been claimed to be young based on near-IR spectra which 
indicate low surface gravity. Examples include the L dwarf 2MASS J01415823-4633574 (Kirkpatrick 
et al. 2010; Patience et al. 2012) and the nearby, 7.5 pc, L5 brown dwarf 2MASS J035523+113337 
(IFaherty et al.ll2012l ). If these 2MASS objects are indeed that young, then their masses might be 
as low as 10 Mj up . As noted above, if WISE 1828+2650 were as young as 0.1 Gyr, its mass would 
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be only 1-2 Mj up . A young, brown dwarf of this age is certainly a possibility, albeit one with a low 
probability (~1%) based on the statistics of nearby low mass objects, i.e. 3 young M stars plus 1 
young L dwarf out of 250 low mass objects. 

The simplest metric for assessing the age of WISE 1829+2650 is the directly measured tangen- 
tial velocity which comes from the proper motion and distance: vtn .n = 4.74/i/II = 51 ± 5 km s _1 
where n is the total proper motion and II the parallax (jSmart II1977I ). This value falls in the middle 
of the distrib ution of tangential v elocities measured for nearby (< 20 pc) L and T brown dwarfs 
(see Figure 7, iFaherty et al.l (120091 )). These velocities are consistent with models of disk population 
with ages of 3-8 Gyr. This range narrows to 2-4 Gyr if the highest velocity objects (v ta n > 100 km 
s _1 are excluded. Within the age range and using fits to the COND models described above, we 
estimate that the mass of WISE 1828+2650 is - 3-6 M Jup . 

A number of authors have used kinematic properties to associate brown dw arfs with nearby 
young clusters, but the results to date are only suggestive (IFaherty et al.ll2012l ). For example, 
might WISE 1828+2650 be a young object could be asso ciated with one of the nearby young 
moving groups, e.g. AB Dor, /3 Pic or Tucana/Horlorgium (IZuckerman fc Songll2004l ; iTorres et al 
20081 ) or the Argus/IC 2391 group as in the case of AP Col? These groups are located at distances 
of 20-100 pc with ages ranging from 8 to 50 MYr. We evaluated possible galactic velocity vectors for 
WISE 1828+2650, (f7, V, W), and their associated uncertainties cruy,w, from its position, parallax, 
and proper motion vector determined here for a range of the unknown V rad from -100 to 100 km 
s _1 . We calculated the difference between cluster {U c u V c i, W c i) values and possible values for WISE 
1828 as a function of V rad \ 



Auvw(v rad ) = V(u W is E (v rad ) - u cl y + (v WISE (v rad ) - v cl y + (v WISE (v rad ) - v cl y (3) 



For V rad between -100 to 100 km s 1 there is no goo d match in (U, V, W) sp a ce for any of 



the six moving groups with ages from 8 to 50 Myr (Table IZuckerman et al. I (120041 ) ; iRiedel et al 
(1201 ll )). The smallest separation is AUVW ~ 50 km s _1 . While it is impossible to rule out the 
possibility that WISE 1828+2650 is a young straggler ejected from one of these groups, it is more 
likely that WISE 1828+2650 is an older, higher mass object. 



5.3. Brown Dwarf or Free Floating Planet? 

The suggestion that a population of free floating planets may be responsib le for short dura - 



tion microlensing events, t ~ y/ 'M /Mj up day ~ 0.3 — 3 day for M=0.1-10 Mj up (jSumi et al.l 120111 ) 
immediately l eads to the c onjec ture that objects like WISE 1828+2650 might be members of such 
a population. ISumi et al.l ()201ll ) calculate that fitting the distribution of short duration microlens- 
ing events requires a population of ^1-2 Mj up objects roughly 1.3-1.9x as numerous as all main 
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sequence stars (0.08-1 M ). Quanz et al. fcoid ) argue that imaging surveys may be starting to 
reveal a significant number of bound planets located at 10s of AU from low mass stars which would 
decrease the need for a large population of free-floating, planetary mass objects. 



The RECONS database fjHenry et al. 112000 ) lists 300 main sequence stars (M-G) within 10 pc, 
or 0.072 object pc -3 . Estimates of the entire population of Y dwarfs down to the coolest objects 
are premature but we note that WISE 1828+2650's distance implies a minimum local density of 
1.7 ± 0.5 x 10 -4 pc -3 . For WISE 1828+2650 to represent the tip of the microlensing iceberg, the 
population of similar objects would h ave to rise very sha rply below the few Mj up mass of this object 
Between 0.1 and 10 Mj up where the ISumi et al.1 (120111 ) data are sensitive (Figure 9 of ISumi et al 



fcoill )). the mass distribution function, dN/DM oc M~ a , would have to have to be quite steep. 
Estimating the mass of WISE1828+2650 using a uniform distribution between 3-6 Mj up yields 
a = 2.2 ± 0.2 to produce enough extremely low mass dwarfs (0.1-10 Mj up ) equal in number to the 
main sequence stars within 10 pc. Thi s is dr amatically steeper than the flat Brown Dwarf index, 
a ~ (Figure 14 in iKirkpatrick et al.1 (120121 )), and steeper than the mass distribution o f planets 
determined by radial velocity studies, a = 1.5 (jCumming et al.ll2008l ; iHoward et al.ll2010l ). 



To date, fourteen "f ree-floating" Y dwarfs have been identified 



20121: rTinnev et al.ll2012r ). Several other WISE-selected Y dwarf candidates, such as those presented 
in lGriffith et al.l ~ 2012 ). remain but lack spectroscopic confirmation. All are on active parallax and 
spectroscopic classification campaigns, so it will soon be known if other WISE 1828+2650 and 
WD 806-661 B analogs exist. Even if, for example, the local density of objects as cool as WISE 
1828+2650 eventually proves to be 5 times higher than estimated solely on the basis of this one 
object, then the power law index required to account for the microlensing results would still be 
> 1.5. In the absence of convincing evidence for an upturn in the mass distribution, the most 
conservative interpretation is that WISE 1828+2650 sits at the low mass end of the brown dwarf 
distribution rather than at the start of a much larger population of free floating "planets". 

A less conservative interpretation of these results is also possible. As Luhman et al. (2011) 
themselves suggest WD 0806-661 B could have formed as a ^6-9 Mj up object in a protoplanetary 
disk around its ^2 M primary, WD 0806-661 A. Only later, after the primary shed mass during its 
evolution off the main sequence, would WD 0806-661 B have migrated into a larger orbit, currently 
of radius ^2500 AU. The formation mechanism between the WISE Y0 and Yl discoveries, which 
presumably formed via cloud fragmentation as do normal stars, and WB 0806-661 B would thus 
be different. Compositional differences should be seen between the two types of objects, given 
that chemical differentiation is expected in a proto-planetary disk and not i n a star formatio n 
cloud. Jupiter's enhanced metallicity is an obvious example of this difference (IWong et al.l 12004 ). 
Objects formed via co re accretion may also have a pronounced core of heavy elements as inferred 
for transiting planets (IRogers fc Seagerll2010l ) which might manifest itself in significant ways. Such 
differences between these two classes could possibly account for the discrepancies noted in the 
absolute magnitude versus spectral type diagrams of Figure 
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Taking this one step further, the similarity between WD 0806-661 B and WISE 1828+2650 may 
be a direct reflection of the fact that both objects were birthed via the same process. In this case, 
WISE 1828+2650 would be a "free-floating" planet that, given our mass esti mate of ~3- 6 Mj 7/ ,„ , 
would fall at the high-mass end of the rogue planet population proposed by ISumi et al.l (120111 ). 
In this case, the poor agreement with theoretical spectra may be a consequence of a poor match 
between the true elemental abundances or interior structure of WISE 1828+2650 and those in the 



model assumptions. If the space density of truly "free-floating" planets measured by ISumi et al. 



(1201 ll ) is correct, then there should be colder examples (<300K) at much smaller distances (<<10 
pc) than WISE 1828+2650 itself. Some of these may yet be discoverable at lower signal-to-noise 
ratios in the WISE data. 



6. Conclusions 

We have used images from four telescopes to determine the parallax and distance to the >Y2 
dwarf WISE 1828+2650. At 11.2^'q pc, the object is among the coolest, lowest luminosity objects 
yet discovered in the solar neighborhood. The high tangential motion of the object suggests an 
age of 2-4 Gyr with a mass of 3-6 Mj up . However, the absolute magnitude of WISE 1828+2650 is 
brighter by several magnitudes (depending on wavelength) than extrapolated from other Y dwarfs 
making the true nature of this source somewhat of a mystery. WISE 1828+2650 is similar in color 
and absolute magnitude to the cool object orbiting the nearby white dwarf WD 0806-661 B. The 
exact nature and evolutionary state, including its mass and age, of WISE 1828+2650 will require 
further observation and theoretical investigation. 

The proximity of WISE 1828+2650 and other Y dwarfs will allow more precise determinations 
of their motions in the coming years as well as a search for any astrometric companions that 
might fall below detectability via direct imaging. The launch of JWST will allow measurement 
of the spectrum of WISE 1828+2650 and other objects with similarly low T e ff from the deep red 
through to the mid-infrared, resulting in a much more detailed physical characterization and a 
better assessment of their ages and masses. Y dwarfs like WISE 1828+2650 and the companion 
WD 0806-661 B will become touchstones for spectroscopic modeling with direct applicability to the 
atm ospheres of bona- fide planets now being measured with transit spectroscopy and direct imaging, 



e.g. iBarman et al.1 (120111 ). The differences expected between bound objects ("planets") and free- 
floating objects (likely brown dwarfs, but possibly free floating "planets") should produce distinct 
spectroscopic markers distinguishable in the high SNR, high spectral resolution measurements 
possible with JWST. 



The research described in this publication was carried out in part at the Jet Propulsion Lab- 
oratory, California Institute of Technology, under a contract with the National Aeronautics and 
Space Administration. This publication makes use of data products from the Wide-field Infrared 
Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet 



-16- 



Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and 
Space Administration. This research has made use of the NASA/IPAC Infrared Science Archive 
(IRS A), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Administration. This work is based in part 
on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract with NASA. Support for this work 
was provided by NASA through an award issued to program 70062 and 80109 by JPL/Caltech. 
This work is also based in part on observations made with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations 
are associated with program 12330. Support for program #12330 was provided by NASA through 
a grant from the Space Telescope Science Institute. Some data presented herein were obtained 
at the W. M. Keck Observatory from telescope time allocated to the National Aeronautics and 
Space Administration through the agency's scientific partnership with the California Institute of 
Technology and the University of California. The Observatory was made possible by the generous 
financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge 
the very significant cultural role and reverence that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community. We are most fortunate to have the opportunity to 
conduct observations from this mountain. In addition we wish to acknowledge the generosity of 
Jay Anderson in sending HST data for the clusters M15 and M92 and useful discussions with Andy 
Gould. The RECONS database of nearby stars remains an invaluable resource. We used the IRSA 
archive at IPAC to access the 2MASS catalogs. 



7. Appendix I. Differential Refraction and Distortion Mapping 



The positions of sources in the Keck images were corrected for the effects of differential refrac- 
tion. The error introduced by coadding before source extraction is small compared to other sources 
of astrometric error. We applied to each extracted position a correction A a j = A a ^(HAo, So, Ao) — 
A a j(HA + AHA, So + AS, Ao) relative to the center of the field where HAq and So are the hour 
angle and declination at the center of the field with the HA specified at the middle of the roughly 
hour-long observational sequence. Ao is the central wavelength of the H-band filter. Figure [9] shows 
2 nd order corrections to the differential refraction as a function of Hour Angle during the duration 
of one observation and for a change in effective wavelength across the H band. The curves show the 
changes in the differential refraction from the nominal 1 st order correction at one corner of the Keck 
field of view (^20v / 2 // , center to corner). At the level of precision required for these measurements, 
these 2 nd order differential effects can safely be ignored over the small Keck field. 

To assess the effects of optical distortion on NIRC2 images in the wide field mode (40 mas 
pixels) we followed the proc edure outlined in | Yelda et al.l (120101 ) using HST observations of the 
cluster M15. The HST data (jAnderson et al.ll2008l ) were corrected for the effects of non-orthogonal 
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axes (coordinate skew). Keck observations of M15 were made in the H band using NIRC2 with a 
natural guide star Adaptive Optics at 2 epochs (Oct 2011 and July 2012). In July 2012 observations 
were made at 3 different orientations of the NIRC2 camera (0°,±45°). At each rotation angle and 
epoch the telescope was stepped through a 9 position pattern across the cluster for a total of 36 
separate observations of some portion of the cluster with each image containing between 50-250 
matches between the Keck and HST source lists. The Keck data were extracted using IDL aper 
routine with a 5 pixel radius. After correction for the effects of differential refraction, lists of sources 
extracted from each image were rotated and shifted into registration with the HST stars. In total the 
HST stars probed 5631 unique (x, y) positions in the NIRC2 focal plane (Figure [TU1). Nominal (a, 5) 
positions were obtained for the Keck sources using plate scale and rotational information derived 
from the image registration. Noise-weighted differences between the Keck and HST positions, 
(Aq, AS), were fitted to a fourth order polynomial as a function of focal plane coordinates (x,y) 
using the Mathematica routine Linear ModelFit. The derived coefficients for transformation from 
native, distorted (x, y) pixel coordinates into un-distorted, (x r , y') pixel coordinates are given in 
Eqns. 4-5 and Table [91 The coordinate system is centered at the center of the NIRC2 array and 
extends from ±512. 



x' =a + a\x + a 2 y + a 3 x 2 + a^xy + a 5 y 2 + a 6 x 3 



- a 7 x 2 y 



agxy + agy + a\0x + a\\x y + a\2x y + a\3xy + y 



y' =6o + bix + b 2 y + b 3 x 2 + b^xy + b 5 y 2 + b 6 x 3 + b 7 x 2 y 
b s xy 2 + b 9 y 3 + &i0x 4 + hlx 3 y + b x 1x 2 y 2 + brfxy 3 - 



(4) 
(5) 



Standard equations then project the (x 1 \y r ) pixel coordinates onto the sky accounting for 
plate scale and rotation angle. O ur H-band plate scale of 39.76 5=b0.004 mas/pixel is consistent 
with the 39.82±0.25 mas/pix from lMetchev fc Hillenbrand! (120041 ) from comparison with a known 
binary system and somew hat smaller than the 39.90±0.08 and 39.86±0.02 K s plate scales derived 



separately in x- and y- by 



Pravdo et al 



from comparison with 17 HST stars. 



Examination of the residuals led to the rejection of <2% of the data, i.e. 100 points that 
deviated significantly in either amplitude or orientation from their neighbors. The first plot of the 
Figure [101 triptych shows Keck-HST positions differences as vectors. The central image shows a 
vector representation of the 4 th order polynomial fitted to the distortion map. The third image 
shows the residuals in the Keck-HST differences after subtraction of the polynomial fit. The 
mean value of the Keck-HST residuals after the model subtraction in focal plane coordinates is 
(x,y) = (0.2 ± 10, 1.0 ± 13) mas. The distortion map was applied to the individual Keck images 
using an IDL script for "de-warping" developed by the Keck Observator}|§. 



http: / / www2.keck.hawaii.edu/inst /nirc2 / forReDoc / post_observing/ dewarp / 
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Table 1. Observing Log 



Telescope 


Epoch 


Modified Julian Date 


Comment 


WISE 


2010-03-30 


55285.00 




Keck/NIRC2-LGS 


2010-07-01 


55378.44 


PI, Wright 


Spitzer/IRAC 


2010-07-10 


55387.34 


PI, Mainzer, PropID=9/551 


WISE 


2010-09-28 


55467.00 




Spitzer/IRAC 


2010-12-04 


55534.27 


PI, Kirkpatrick, PropID=10/70062 


HST/WFC3-IR 


2011-05-09 


55690.89 


PI, Kirkpatrick,PropID=12330 


Keck/NIRC2-LGS 


2011-10-16 


55850.21 


PI, Beichman 


Spitzer/IRAC 


2011-11-19 


55894.04 


PI, Kirkpatrick,PropID= 10/80109 


Spitzer/IRAC 


2012-05-25 


56072.25 


PI, Kirkpatrick,PropID= 10/80109 


Keck/NIRC2-LGS 


2012-07-09 


56117.32 


PI, Beichman 


Keck/NIRC2-LGS 


2012-10-07 


56207.22 


PI, Beichman 
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Table 2. HST Sources in Other Datasets 



RA 


DEC 


F140W 


a 


Center 


Posn. 


Dataset 


(J2000) 


(J2000) 


(mag) 


(mag) 


Dist.(") 


Err. (mas) 




277.10529566 


26.85705074 


15.66 


0.001 


84.9 


42.3 


Spitzer, 2MASS 


277.10824423 


0/° OTA0001 O 

26.85082312 


in r r> 

19.56 


n non 

0.020 


f* P* <"7 

65.7 


o o o 

32.8 


Spitzer 


277.11519519 


26.83369321 


1 >~7 OT 

17.27 


r\ nno 

0.002 


r o f 

52.5 


26.2 


Spitzer 


277.11647015 


26.84300085 


16.57 


r\ nno 

0.002 


O A O 

34.2 


17.0 


Spitzer, 21V1ASS 


277.11790606 


O/^ or'inor'p'O 

26.85103552 


19.63 


0.009 


39.6 


19.7 


Spitzer 


277.12026174 


i a /■ i o o o o n /i p" z' 1 

26.82339456 


17.37 


0.002 


76.3 


38.0 


Spitzer 


277.12155151 


26.84341695 


18.21 


0.003 


17.8 


8.9 


Keck 


277.12463494 


26.84734988 


19.24 


r\ r\r\ f 

0.005 


1 P* o 

15.3 


7.7 


Keck 


277.12479343 


26.84571497 


21.22 


0.018 


10.3 


5.6 


Keck 


O^T^T lO/in/l/^OO 

277.12494638 


26.84401162 


20.30 


0.010 


7.0 


O i~7 

3.7 


Keck 


on , '7 io/^/°0/ion 

277.12662489 


26.84369877 


21.76 


0.025 


1.5 


3.3 


Keck 


277.12761710 


26.86279933 


17.29 


0.002 


68.8 


34.2 


Spitzer 


ott i OTnnoo /I 

277.12799324 


26.84048725 


19.13 


0.005 


11.9 


6.0 


Keck 


277.12816179 


26.84302297 


20.64 


0.013 


A O 

4.2 


2.7 


Spitzer, Keck 


c%tt i o o /° p" /° p" 

277.12865656 


o/° o o nn o t p" p" 

26.83903755 


18.24 


0.003 


17.5 


8.7 


Keck 


CiTT 1 00^7^700^7 

277.12877327 


26.85757557 


15.28 


0.001 


50.2 


25.0 


o • i oi\ /r A o o 

Spitzer, 2MASS 


277.12906930 


o/° o At nn/°o /) 

26.84100634 


in f o 

19.52 


r\ r\r\r' 

0.006 


ii/^ 
11.6 


f o 

5.8 


Keck 


277.12971270 


o/° o r r o /i /° ^ i 

26.85534671 


it r n 

17.50 


r\ nno 

0.002 


42.7 


1 o 

21.3 


Spitzer 


277.13008161 


0/° O A AA r 7 A r r 

26.84907455 


on f f 

20.55 


0.011 


O 1 o 

21.6 


1A n 

10.8 


Spitzer, Keck 


277.13061562 


26.84941597 


20.97 


0.016 


23.5 


11.9 


Spitzer, Keck 


277.13116667 


26.84513525 


20.09 


0.008 


i/ii 
14.1 


^7 1 

7.1 


Keck 


277.13117765 


26.85026363 


19.96 


0.008 


27.0 


13.5 


Spitzer, Keck 


277.13143360 


O/^ O/IOOO/^OO 

26.84232633 


22.04 


0.034 


14.8 


O f 

8.6 


Keck 


c\i~ri~7 i 01 p"oono 

277.13152398 


26.84961178 


20.11 


0.008 


25.6 


12.8 


Spitzer, Keck 


277.13173964 


OP"/!Oni^7P"' 

26.85489175 


17.43 


0.002 


42.9 


21.4 


Spitzer 


C\<~7<~7 1 OOI n n /I 

277.13310074 


26.86055572 


16.21 


0.001 


63.7 


31.7 


Spitzer, 2MASS 


277.13366038 


26.82526749 


i n oo 

16.28 


r\ nno 

0.002 


69.6 


O A >"7 

34.7 


o • X oi\ /r A o o 

Spitzer, 21V1ASS 


277.13461287 


0/° OO/IO/IOOO 

26.83484333 


i o or 

18.35 


r\ AA A 

0.004 


40.0 


inn 

19.9 


Spitzer 


flTT 1 O^/^OO/IO 

277.13763243 


0/° 001/°10/10 

26.83161343 


i /"» on 

19.20 


r\ nn ^7 

0.007 


f f i 
55.1 


27.5 


O ' 4- 

Spitzer 


277.14074059 


26.83242086 


16.24 


0.001 


59.8 


29.8 


Spitzer,2MASS 


277.14122122 


26.84312777 


17.25 


0.002 


45.4 


22.6 


Spitzer 


277.14136220 


26.84167755 


17.19 


0.002 


46.4 


23.1 


Spitzer,2MASS 


277.14168068 


26.83999147 


15.44 


0.001 


48.7 


24.3 


Spitzer,2MASS 


277.14393913 


26.84815732 


16.63 


0.001 


56.4 


28.1 


Spitzer,2MASS 


277.14408060 


26.83259772 


14.24 


0.001 


67.6 


33.7 


Spitzer,2MASS 


277.14510022 


26.83998378 


17.04 


0.002 


59.4 


29.6 


Spitzer 


277.14677186 


26.83408275 


14.78 


0.001 


72.1 


35.9 


Spitzer,2MASS 


277.14732238 


26.83496895 


17.38 


0.002 


72.2 


35.9 


Spitzer 
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Table 3. Photometry for WISE 1828+2650 



Telescope 


Date,UT 


Filter 


A(/im) 


Maj 


y. (mag) 


SNR 


F^uJv) 1 


Keck/NIRC2 


2010-07-01 


J 3 


1.25 


23 


57±0.35 


3 


0.59± 0.21 


HST/WFC3 


2011-05-09 


F140W 


1.39 


23 


15±0.09 


11 


0.75± 0.07 


Keck/NIRC2 


Average 


H 


1.65 


22.45±0.08 2 


17 


1.07± 0.09 


Keck/NIRC2 


2010-07-01 


H 


1.65 


22 


52±0.23 


5 


1.34± 0.25 


Keck/NIRC2 


2011-10-16 


H 


1.65 


22 


42±0.14 


11 


1.24± 0.13 


Keck/NIRC2 


2012-07-09 


H 3 


1.65 


22 


46±0.12 


11 


1.40± 0.12 


WISE 


2010-03-30 


Wl 4 


3.35 




>18.5 




<12 


Spitzer/IRAC 5 


2010-07-10 


Chi 


3.55 


17 


17±0.03 


55 


38.1± 1.1 


Spitzer/IRAC 


Average 


Ch2 


4.49 


14 


66±0.03 


75 


247± 9 


Spitzer/IRAC 


2010-07-10 


Ch2 


4.49 


14 


61±0.05 


60 


256± 13 


Spitzer/IRAC 


2010-12-04 


Ch2 


4.49 


14 


60±0.05 


40 


261± 13 


Spitzer/IRAC 


2011-11-29 


Ch2 


4.49 


14 


70±0.05 


40 


238± 12 


Spitzer/IRAC 


2012-05-25 


Ch2 


4.49 


14 


72±0.05 


50 


232± 12 


WISE 


2010-03-30 


W2 


4.6 


14 


39±0.06 


20 


299± 18 


WISE 


2010-03-30 


W3 4 


11.56 




>12.5 




<290 


WISE 


2010-03-30 


W4 4 


22.08 




>8.8 




<2501 



Note. — 1 Flux densities have not been color-corrected. 2 Uncertainty weighted av erage, including 
5% calibration uncertainty. 3 Keck NIRC2/LGS result from lKirkpatrick et al.l (1201 ll ). 4 Upper limits 
are 95% confidence. 5 Includes 3% IRAC calibration uncertainty. 
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Table 4. Registration of Images 



Tel-l/Tel-2 


Date(UT) 


Number of 


Aa 1 


AS 1 


ff meoti 


Rot. Unc. 3 






Stars 


(mas) 


(mas) 


(mas) 


(deg) 


#ST/2MASS 4 


2011-05-09 


11 


10 


-2 


±4 


0.028 


Keck/ 'HST 


2010-07-01 


11 


0.0 


0.0 


±3 


0.003 


Spitzer/HST 


2010-07-10 


22 


2 


-11 


±11 


0.010 


Spitzer/HST 


2010-12-04 


21 


17 


11 


±16 


0.017 


Keck/HST 


2011-10-16 


21 


1.7 


0.0 


±3 


0.003 


Spitzer/HST 


2011-11-29 


16 


-22 


17 


±18 


0.019 


Spitzer/HST 


2012-05-25 


15 


5 


-16 


±13 


0.015 


Keck/HST 


2012-07-09 


15 


-1.8 


2.4 


±3 


0.003 


Keck/HST 


2012-10-07 


20 


-0.2 


-0.1 


±12 


0.009 



Note. — 1) Difference between two reference frames; 2) Standard deviation of the mean, 
(a/y/N),in the offset in RA, DEC between stars observed with both telescopes; 3) Uncertainty 
in image rotation; 4) after removal of (Aa , AS) (initial HST-2MASS)=(-0.50±0.02,0.15±0.05) // off- 
set. 
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Table 5. Astrometric Data for WISE 1828+2650 



Telescope 


Date(UT) 


MJD 1 RA (sec) 2 


DEC(") 2 


UncertaintyC") 3 




WISE 4 


2010-Mar-30 


55284.50 


31.07805 


37.9040 


0.210 




Keckl 


2010-Jul-01 


55378.44 


31.10325 


37.9443 


0.029 




Spitzerl 


2010-Jul-10 


55387.34 


31.10961 


37.9780 


0.058 




WISE 4 


2010-Sep-28 


55467.00 


31.09463 


37.8860 


0.210 




Spitzer2 


2010-Dec-04 


55534.27 


31.12461 


37.9020 


0.062 




Hubble 


2011-May-09 


55690.89 


31.17538 


38.0725 


0.030 




Keck2 


2011-Oct-16 


55850.21 


31.19329 


38.0785 


0.025 




Spitzer3 


2011-Nov-29 


55894.04 


31.20339 


38.1592 


0.062 




Spitzer4 


2012-May-25 


56072.25 


31.24935 


38.2111 


0.061 




Keck3 


2012-Jul-08 


56117.32 


31.25462 


38.3020 


0.025 




Keck4 


2012-Oct-07 


56207.22 


31.25951 


38.1087 


0.040 




Note. - 


1) MJD=Modified Julian Date; 2) Relative to 18 h 


28 m and +26°50 m (J2000.0): 


; 3)Ex- 


eludes ^80 mas uncertainty 


in absolute 2MASS reference frame. 4) Not included in fitting 


due to 


large uncertainties. 












Table 6. Astrometric Data for WISE 1828+2650 With Respect to Source A 




Telescope 


MJD 1 Aa 


(BD-SrcA,") 


A5(BD-SrcA,") Uncertainty^) 2 




Keckl 


55378.44 


4.647 




3.042 


0.032 




Spitzerl 


55387.34 


4.609 




3.179 


0.121 




Spitzer2 


55534.27 


4.896 




3.009 


0.122 




HST 


55690.89 


5.570 




3.192 


0.018 




Keck2 


55850.21 


5.853 




3.175 


0.015 




Spitzer3 


55894.04 


5.884 




3.264 


0.134 




Spitzer4 


56072.25 


6.535 




3.301 


0.134 




Keck3 


56117.32 


6.615 




3.374 


0.016 




Keck4 


56207.22 


6.695 




3.303 


0.042 





Note. — 1) MJD=Modified Julian Date; 2) Includes uncertainties due to extraction uncertainties, 
field rotation, residual Keck distortion. 
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Table 7. Parallax And Proper Motion 



Parameter Method-1 1 Method-2 2 Method-(l+2)( 2 ) 



("/yr) 


0.969±0.017 


0.936±0.015 


0.954±0.011 


us ("M 


0.157±0.019 


0.152±0.017 


0.153±0.0125 


Parallax, 7r,(") 


0.103±0.016 


0.079±0.012 


0.090±0.0095 


Distance (pc) 




12.6±?;| 


11 2 +L3 
iJ -- z -i.o 


"5/2000.0 


30.344 s ±0.077 s 


30.375 s ±0.012 s 


30.359s ± 0.009s 


ri 3 

°J2000.0 


36.120" ±0.18" 


36.200" ± 0.21" 


36.278" ±0.147" 


X 2 (with # deg. of freedom) 


18.6 (13) 


12.7 (13) 


36.3 (31) 


Adopted Distance 4 (pc) 


11 2 +L3 






Tangential Velocity 4 , vt an (km s _1 ) 


4 r+15 







Note. — 1) Based on "absolute coordinates"; 2) Based on offsets from Source A; 3) With respect 
to 18^28 m and 26°50 / . The random uncertainty does not include ^80 mas uncertainty in absolute 
2MASS reference frame. 4) From the combination of Method-1 and Method-2. 



Table 8. Cluster Kinematics and WISE 1828+2650 



Name 1 


Age 


IJ(2) 


V 


W 


Min AUVW { V 




(Myr) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


Argus/IC2391 


40 


-22 


-14.4 


-5 


46 


TW Hydrae 


8 


-11 


-18 


-5 


50 


Tucana/Horologium 


30 


-11 


-21 





55 


f3 Pictoris 


12 


-11 


-16 


-9 


48 


AB Doradus 


50 


-8 


-27 


-14 


49 


r] Chamaeleontis 


8 


-12 


-19 


-10 


46 


Columba 


30 


-13 


-22 


-6 


50 



Note. — 1 Cluster data from IZuckerman et al. I (120041 ) : iTorres et all (120081 ) : 2 Cluster velocities 
are calculated such that U is positive toward the Galactic Center; 3 Smallest difference between 
([/, V, W) for WISE 1828 and a cluster within a range of unknown V ra d between -100 and 100 km 
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Table 9. H Band Distortion Coefficients 



Parameter X-Estimate X-Uncertainty Y-Estimate Y-Uncertainty 



const. 


-0.306229 


6.16x10" 


-4 


-0.700962 


6.16x10" 


-4 


X 


0.00369728 


2.44x10" 


-6 


-0.00134633 


2.44x10" 


-6 


y 


-0.0021500 


2.39x10" 


-6 


-0.00333556 


2.39x10" 


-6 


x 2 


7.931 xlO" 7 


1.17x10" 


-8 


6.116) 


<10" 6 


1.17x10" 


-8 


x y 


-1.112X10" 6 


1.02x10" 


-8 


2.258) 


<10" 7 


1.02x10" 


-8 


y 2 


1.0698xl0" 5 


1.15x10" 


-8 


5.081) 


<10" 6 


1.15x10" 


-8 


X 3 


-5.857xl0" 9 


1.47x10" 


-11 


-1.957X10" 10 


1.47x10" 


-11 


x 2 y 


1.015xl0" 9 


1.39x10" 


-11 


2.677) 


<10"9 


1.39x10" 


-11 


xy 2 


-3.869xl0- 9 


1.35x10" 


-11 


1.060) 


<10"9 


1.35x10" 


-11 


y 3 


6.10X10" 10 


1.36x10" 


-11 


5.217) 


<10"9 


1.36x10" 


-11 


x 4 


-8.486 xlO" 12 


5.21x10" 


-14 


9.835x 


10" 13 


5.21x10" 


-14 


x 3 y 


9.712xl0" 13 


4.90x10" 


-14 


-3.382) 


<10" 12 


4.90x10" 


-14 


x 2 y 2 


-5.878xl0" 12 


4.73x10" 


-14 


2.080x 


lO" 12 


4.73x10" 


-14 


xy 3 


2.695xl0" 12 


4.78x10" 


-14 


-4.396) 


<10" 12 


4.78x10" 


-14 


y 4 


-2.244xl0" 12 


4.88x10" 


-14 


-2.219) 


<10" 12 


4.88x10" 


-14 



-28- 




^ - 














QD 










31.: 

(A 
CN 
lO 


Is 31. 


s 18h28m: 


fl%s 30.1 


)s 30. 


s 


E 
8 

■o 

(0 














32$s 



31.9s 



18h28m<0.9s 



29.9s 




1 



28.9s 



Fig. 1 — An HST WFC3/F140W image of WISE 1828+2650 with the reference stars used for the 
co-registration of the Spitzer, HST, and Keck fields circled in red; a few of the 2MASS reference 
stars fall outside of the image shown here to preserve an adequate scale in the vicinity of the WISE 
source. The absolute positions were derived with respect to 11 2MASS stars spread over the entire 
HST field. The positions of WISE 1828+2650 at different epochs are shown in an inset with HST 
positions in red, Spitzer positions in green, and Keck positions in orange. The alternate astrometric 
reference denoted "Source A" appears to be slightly extended and is located to the southwest of 
WISE 1828+2650. 
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Fig. 2.— A Keck-II H-band image of WISE 1828+2650 obtained with the NIRC2 instrument using 
the laser guide star Adaptive Optics system. The image is 25" across with North up and East to the 
left. The solid line is 2" in length and sits just to the East of WISE1828+2650 which is circled in red. 
Also shown are three of the reference stars used to establish the astrometric reference (Figure [1]). 
The alternate astrometric reference ( "Source A" ) is circled and is located to the southwest of WISE 
1828+2650. 
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Fig. 3.— The fit to the positions of WISE 1828+2650 obtained using Method-1 over 2 years of 
observations. Differences in right ascension (Aa cos (5), upper blue lines and circles) and declination 
(A5, lower red lines and squares) between WISE 1828+2650 relative to its initial WISE position 
are shown. Data points are labeled with their appropriate observatory. The solid lines represent 
the proper motion and parallactic effect as seen from an observatory located on the Earth or in 
Earth-orbit while the dashed lines represent the view from Spitzer's drift-away orbit. 



-31 - 




Fig. 4. — The fit to the positions of WISE 1828+2650 obtained using Method-1 over 2 years of 
observations with the effects of proper motion removed show the fit to the parallactic ellipse. The 
inner and outer dashed ellipses (blue) show the ±1 a parallax values. Individual sightings are 
labeled sequentially for each telescope and have an arrow pointing to their expected location on 
the parallactic ellipse for the time of observation. 
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Fig. 5. — The histogram (left) represents the results of a 10,000 iteration Monte Carlo calculation for 
the combined Method- 1 and Method-2 datasets wherein the distance solutions were generated from 
synthetic data generated with a Gaussian distribution with appropriate mean and uncertainties 
(Tables [Sj). The cumulative distribution function on the right gives derived distances and was used 
to define the uncertainties in the distance. 
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Fig. 6. — The spectral energy distribution for the Y dwarf WISE 1828+2650 inclu des data from 



HST : Keck, WISE and Spitzer (black diamonds). The COND model calculation (IBaraffe et al 



20031 ) for a T e fj = 350 K object is shown in red, bracketed by the predicted emission in the various 



bands for (0.1 Gyr,0.5 M Jup , 240 K) and (10 Gyr,20 M Jupi 389 K) objects with green and blue 
circles, respectively. 



-34- 




Fig. 7. — The quality of fit (x 2 ) between the photometric observations (Keck H-band, HST F140W, 
Spitzer channels 1 and 2) and the COND models is shown as a function of Age (Gyr) and Mass 
(Mj up ). The contours give the logarithm of the x 2 °f the fit for 2 degrees of freedom (4 data points 
minus two derived values). As discussed in the text, the best fit is for a very young, low mass 
object which is not likely on dynamical grounds. For an adopted dynamical age of 2-4 Gyr, the 
best fitting mass is in the range of 3-6 (Mj up ). 
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Fig. 8. — Absolute magnitude versus spectral type diagrams for late-T and Y dwarfs at the H-band 
(left) and WISE W2-ban d (right) for objects with secure parallaxes. The figures are adapted from 



Kirkpatrick et al.1 (120121 ) and include revised parallax me asurements as repor ted in iMarsh et al . 
(20]j). Anew Y dwarf (WISE J163940.83-684738.6) from lTinnev et al.1 (El) has been added to 



the Mw2 vs type plot. Shown by the red, solid point is the location of WISE 1828+2650. Dashed 
red lines mark the absolute magnitude upper limit at H-band and the one-sigma upper and lower 
limits in the W2-band for WD 0806-661B. Because there is no spectrum for WD 0806-661B, the 
spectral type range shown is a n estimate. (For this source, the H magnitude limit is based on 
the J magnitude limit given by iLuhman et al.l (120121 ). assuming the J-H color of 0.7 from WISE 
1828+2650; WISE does not detect this source, so the IRAC channel 2 magnitude, which is nearly 
equivalent to the W2 magnitude, is plotted instead.) Note that at W2 in particular, both WISE 
1828+2650 and WD 0806-661B lie significantly above the trend established for earlier Y dwarfs. 
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Fig. 9. — Effects of 2 order differential refraction as a function of hour angle and central 
wavelength for observations of 1828+2650. 
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Fig. 10. — NIRC2 -LGS observations of the cluster M15 were compared with observations taken 
with HST/WFC3 (j Anderson et allbood ). left) The difference between the HST and Keck position 
as a function of position in the Keck focal plane; center) a 4th order polynomial fit to the distortion 
map; right) residuals remaining after subtracting the model from the data. 



